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Abstract. We report a first in modeling and simulation of the effects of the HIV proteins on the 
(caspase dependent) apoptotic pathway in infected cells. This work is novel and is an extension 
on the recent reports and clarifications on the FAS apoptotic pathway from the literature. We 
have gathered most of the reaction rates and initial conditions from the literature, the rest of the 
constants have been computed by fitting our model to the experimental results reported. Using 
the model obtained we have then run the simulations for the infected memory T cells, called also 
latent T cells, which, at the moment, represent the major obstacle to finding a cure for HIV. We 
can now report that the infected latent T cells have an estimated lifetime of about 42 hours from 
the moment they are re-activated. As far as we know this is the first result of this type obtained 
for the infected memory T cells. 

1 Introduction 

Apoptosis (or programmed cell death) is the process through which a multi-cellular organism 
removes the cells that pose a threat locally (at the tissue level) or globally to the organism. 
During the development of organisms, apoptosis is initiated in order to remove unneeded cells at 
various developmental stages. Apoptosis is a form of cellular suicide, an attempt by an organism 
to cleanly rid itself of cells that are unwanted, improperly formed, injured, infected, etc. It is a 
vital process of multicellular organisms, because of this, the mechanisms of apoptosis are hard- 
coded in all the cells of an organism and it can be executed in several ways; we will focus on 
one particular signaling pathway, which utilizes members of the Caspase protein family. The 
apoptotic pathway has been recently under a high degree of research focus |19|35|36] due to the 
fact that it is estimated that more than 50% of the cancers result due to the inhibition or damage 
in apoptotic signal transduction. Using the model from [19] we have recently simulated the FAS- 
mediated apoptosis in [9] by using a deterministic waiting time algorithm from [26] which was 
then modified and improved into a nondeterministic waiting time algorithm in [21]. The main 
benefit for using the nondeterministic algorithm proposed by us is the fact that it is a discrete 
simulation method (accounting for the discrete nature of the signal transduction pathways) as 
opposed to differential methods and is also fast as opposed to the stochastic simulations such as 
Gillespie's algorithm (see |15|16j ). without losing the nondeterministic behavior. 

We are extending in the current paper the model from [21] to include HIV-1 proteins in the 
nondeterministic discrete simulation of the effects of the HIV on the Fas-induced apoptosis. In 
the following we will give a brief reason for the extension. 

The HIV virus is currently causing a global epidemic according to the report published in 
Geneva on November 21, 2006 by the World Health Organization. The virus infected roughly 
1% of the global population with about 7.3% of the infected people dying in 2006. 
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The HIV-1 virus, which is responsible for 99% of the HIV infections overall, has some re- 
markable features. It infects cells from the immune system of the body: mainly CD4+ T cells, 
but also monocytes and macrophages. The initial infection seems to be through the R5 virus 
type that uses the co-receptor CCR5 for cell entry as reported in [33j. Reports have shown that 
the R5 evolves through mutations into the X4 variant of HIV-1 (using co-receptor CXCR4) in 
the infected individuals late during the course of the infection (see for example |33|11|27] ) . There 
is a significant difference between the R5 and X4 variants of the virus, R5 is inducing apoptosis 
in the bystander T cells while the X4 is responsible for the synctium formation (a multi nucleic 
cell) which proceeds to the killing of even more immune cells. The fact that the X4 is usually not 
present in the first stages of the HIV infection, but is more prevalent in the last stages of AIDS 
suggests that a normal immune system could counteract the X4 strand effectively. So one could 
derive that the immune system is first weakened by the R5 strand through apoptosis before the 
X4 strand can "survive" and collapse it. 

In this paper, we focus on simulating the effects on the Fas-induced apoptosis caused by 
the HIV-1 strand R5 proteins. We report a simulation for the Fas- mediated pathway in HIV-1 
infected cells. There are no other simulations of this type in existence. The HIV uses cellular 
machinery to produce its proteins, in the creation of new virions. Therefore, when a T cell is 
activated by the immune system, it can become infected by HIV, produce HIV proteins (and new 
virus') for a couple of days, then die by apoptosis. Yet, sometimes the CD4+ T cells will become 
inactive after the HIV genome is incorporated into the cellular DNA. This results in a so-called 
latent infection, which has been of major focus in recent research (reviewed in [Uj). In [T7] 
the authors explain that latently infected cells are persistent in patients undergoing active HIV 
treatment (such as HAART). It is clear that a better understanding of the latent cells will assist 
biologists in the development of new treatments. In Section 2 we will highlight the key aspects 
of our Fas-induced apoptosis model, Section 3 will describe the new simulation technique. In 
Section 4 we review the role of HIV proteins in cell death and in Section 5 we present the results 
our from model's simulations and we finish with discussions in Section 6. 

2 Fas-mediated signal cascade 

There are several distinct death receptors, which, when activated, can lead to cellular apoptosis 
through a tightly regulated molecular signaling cascade [3] . We have focused our previous efforts 
on simulating the well- understood Fas- mediated pathway |21| . According to the literature [28] 
and [20] . understanding the complex signaling cascade of Fas-mediated apoptosis can assist in 
the creation of remedies for cancer and autoimmune disorders. In this paper we provide an 
expansion of our previous model, described in [21| . in order to decipher and predict the behavior 
of HIV-infected cells. But first, we discuss the Fas- mediated signaling pathway. 

There are two versions of the Fas-mediated pathway [32] , they are generally referred to as the 
type I and type II pathways. Both pathways begin with the Fas ligand binding to the Fas receptor 
(CD95) on the cell membrane. The ligand forces a conformational change of the Fas receptor, 
which then goes on to recruit the Fas-associated death domain (FADD). Once FADD is bound 
(up to three FADD per Fasl-Fas site), Procaspase 8 and FLIP can be competitively recruited to 
the death domain. FLIP is an inhibitor of the pathway, in that its binding to the death domain 
removes binding sites for Procaspase 8. The recruitment of at least two Procaspase 8 molecules 
can lead to cleavage into active Caspase 8 form. The divergence in the two pathways occurs after 
active Caspase 8 is in the cytoplasm; the end result of both pathways is the cleavage of Procaspase 

3 into active Caspase 3. If the initial concentration of Caspase 8 is large enough, Caspase 3 is 
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directly phosphorylated by Caspase 8 (type I pathway). Otherwise, Caspase 8 truncates Bid, 
which binds with two Bax molecules enabling the release of cytochrome c from the mitochondria 
(type II pathway). Bcl-2 can hinder the release of cytochrome c (thus, inhibiting apoptosis) by 
binding to Bid, tBid, Bax or the tBid:Bax molecular complexes. Once released, cytochrome c 
binds to Apaf and ATP before recruiting and activating Caspase 9. The active Caspase 9 can 
then go on to activate Procaspase 3. XIAP can prohibit apoptosis by binding to Caspase 9, thus 
preventing the phosphorylation of Caspase 3. 

Though the molecules of the Fas-mediated pathways (type I and type II) have been analyzed 
experimentally (pro- and anti-apoptotic), we believe the only way to understand the dynamics of 
the complex apoptotic system is through computer simulations. In [19] and [21], three different 
modeling methods are employed that yield results similar to experimental observations. We use 
our ideas from [21] and modify the model to incorporate the functions of HIV proteins. 

3 The membrane system model 

There are three schools of thought on modeling molecular interactions. Microscopic chemistry 
- that is, the consideration of molecular dynamics for every single molecule in the system (po- 
sition, momenta of atoms, etc.) - is computationally intractable; the time scale involved in the 
simulations, uncertainty about cellular components, and the number of atoms involved are fac- 
tors that make this approach impossible to implement in a computational system. Differential 
equations are the macroscopic approach. This approach has been applied to the Fas pathway 
previously in [19]. 

We believe that differential equations are not the best approach for modeling signal cascades 
involving low numbers of molecules. We discussed this in general for the Fas-mediated pathway 
in [21] and now we expand our ideas into the realm of HIV proteins, where proteins like Vpr 
are at a level of only ~700 molecules upon viral infection [6j. Our approach, referred to as 
mesoscopic chemistry, is similar to the microscopic approach. However, we ignore insignificant 
molecules, like water and non-regulated parts of the cellular machinery, and we do not directly 
model the position and momenta of the molecules, but we model molecular interactions using 
statistical equations governing how often and which molecules interact. Our approach allows finer 
modeling of a molecular system than the macroscopic approach, and yet it is more tractable 
than the microscopic approach. 

The membrane system paradigm is based on the compartmental structure found in cells. 
Each compartment has different rules and objects, and the system evolves as the rules are 
carried out based on the objects in the compartments. As we are modeling molecules in a cell, 
it is a natural approach to use a membrane system. 

We provide now all of the different types of rules (reactions) our membrane system can 
simulate: 

Transformation, complex formation and dissociation rules: 



(5) [ a ]i — ► [ b, c, d ]i where a,b,c,d € U and I G L. 

These rules are used to specify chemical reactions taking place inside a compartment of type 
/ € L, more specifically they represent the transformation of a into b; the formation of a 
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complex c from the interaction of a and b; and the dissociation of a complex a into b and c. 
The rules of type (4) model the formation of a complex d from the interaction of 3 molecules: 

a, b and c and the rules of type 5 model the dissociation of the complex a into 3 molecules 

b, c, and d. 
Diffusing in and out: 

(6) [a],-a[], 

(7) a [ ]i — » [ a ]i where a € £ and / € L. 

When chemical substances move or diffuse freely from one compartment to another we use 
these types of rules, where a moves from or to a compartment of type /. 
Binding and unbinding rules: 

(8) o[&],-[c], 

(9) [ a ]i — > b [ c ]i where a,b,c G U and I € L. 

Using rules of the first type we can specify reactions consisting in the binding of a ligand 
swimming in one compartment to a receptor placed on the membrane surface of another 
compartment. The reverse reaction, unbinding of substance from a receptor, can be described 
as well using the second rule. 
Recruitment and releasing rules: 

(10) a[&]i->c[], 

(11) c [ ]i — > a [ b ]i where a, 6, c G S and I € L. 

With these rules we represent the interaction between two chemicals in different compart- 
ments whereby one of them is recruited from its compartment by a chemical on the other 
compartment, and then the new complex remains in the latter compartment. In a releasing 
rules a complex, c, located in one compartment can dissociate into a and b, remaining a in 
the same compartment as c, and b being released into the other compartment. 

Our interest in the simulator lies in the molecular multiplicity evolution over time. Our rules 
are carried out at varying lengths of time, obeying the Law of Mass Action - that is, the reac- 
tion rate depends proportionally on the multiplicities of the reactants. Therefore, the number of 
reactants for a given reaction logically governs the speed at which that reaction will occur. In 
our model, the time for each reaction to occur is calculated (Waiting Time or WT), and after 
sorting all waiting times, the reaction with the least waiting time is applied once (reactants de- 
creased/products increased by one molecule). The simulation clock is aggregated (r = r+WT), 
and the WT of each the reaction affected by the applied reaction (i.e., number of reactants of 
reaction changed) are recalculated. The process continues until the desired simulation time has 
been reached. 

To calculate the waiting time for each reaction r we have an associated reaction rate con- 
stanJE For each reaction r we pre-compute a kinetic constant const r : const r = V i^i k ^ Ni ^i where 
V is the volume of the system, N is Avogadro's constant (6.0221415 x 10 23 ) and i is the number 
of reactants involved in the reaction. The time needed for the reaction r to finish will then be 
computed in the following way: wt r = const ,1 \ aub \ where |A| and \B\ represent the multiplicities 
of the two reactants A, B of the reaction r. On the other hand, if the reaction simulated is a 
dissociation or other first order reaction involving one reactant, A, then the time is computed 
as wt r = k . Now we shall discuss how the HIV proteins interact with the Fas- mediated 
pathways (type I and type II), and the method with which we set up a simulation of these 
effects. 

1 Some of the reaction rate constants are from [IS] and the HIV protein associate rates are derived from various 
sources 
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4 The roles of HIV proteins in apoptosis 



The onset of AIDS is associated with the loss of immune function, due to the death of CD4 T 
cells (among others). As part of our simulation, we have focused on the mechanisms through 
which an infected cell is kept alive in order to produce new provirus', aiding in the spread of the 
infection. 

The method through which an HIV virus infects a cell is well understood. A spike on the 
virus (the gpl20 molecule) binds to the CD4 receptor of primary T cells, and in conjunction with 
subsequent binding to a coreceptor (CCR5 or CXCR4), a path is opened for the virus to inject 
its contents into the cell [5] and [35]. Reverse transcriptase translates the RNA into DNA and 
the code is implanted into the cells own DNA for future production, once the cell is activated. 
During this time, the immune system fails to detect and destroy the infected cell. HIV proteins 
that were packaged in the virion with the RNA interact with various signaling cascades (e.g., Fas 
pathway) within the cell preparing it for viral reproduction while avoiding annihilation by the 
immune system. It is not until the immune system activates the T cell that viral reproduction 
occurs. 

In this study, we focus on the specific functions of the HIV proteins related to the apoptosis of 
CD4 T cells; reviews are found in [31 j . |29j and [2j. The functions of some of these proteins are still 
up for debate, however, we have pooled the collective knowledge of the biological community in 
order to categorize and model the described functions of various HIV proteins. For an illustration 
of the FAS-mediated apoptosis pathway and the involvement of the HIV proteins, we refer the 
reader to Fig. [TJ 

The HIV protein, Tat, has shown both pro-apoptotic and anti-apoptotic behavior. In [3], the 
authors were able to show that Tat expression in the cell leads to up-regulation of Procaspase 
8 (proapoptotic molecule). Tat has been associated with upregulation of Bcl-2 (antiapoptotic 
molecule) [TJJ. Conversely, [31] argues that Tat downregulates Bcl-2, resulting in increased sen- 
sitivity to the Fas-mediated type II pathway. Tat has also been implicated in the upregulation of 
the Fas ligand in infected cells [4] and [40]. However, the upregulation is assumed to be a mech- 
nism to destroy CTLs (thus, keeping the cell alive) and bystander cells through Fas-mediated 
apoptosis, the increase of Fas ligand does not appear to affect apoptosis of the infected cell |31j . 

Another HIV protein, Vpr, has been implicated in both prohibiting and assisting apoptosis. 
Since Vpr is located in the virions, a low concentration is avaible early in the infected cell life- 
cycle (assumed to be 700 molecules according to [6J). At low multiplicity, Vpr has been shown 
to prohibit apoptosis by up-regulating Bcl-2 and down-regulating Bax (proapoptotic molecule) 

M- 

HIV protease, which is responsible for cleaving the proteins created from the viral genetic 
code, has also been shown to affect the signaling cascade in apoptosis. The HIV protease has 
been shown to cleave Bcl-2 into a deactive state |34| . while it is also responsible for cleaving 
Procaspase 8 into active Caspase 8 form [24], both events lead to an increased susceptibility for 
Fas- mediated apoptosis. 

The HIV RNA travels to the nucleus and becomes incorporated into the cellular genetic code. 
While the T cell is activated by the immune system, the HIV proteins are manufactured through 
the cellular machinery (transcription into HIV RNA and translation through the ribosomes in 
the cytoplasm). The simulator makes stochastic choices on where to send the Tat, Vpr, and HIV 
Protease (cytoplasm, mitochondria, or extracellular). Using the information on the three HIV 
proteins, we were able to create 11 new rules for our simulator. Our results are discussed below. 
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Fig. 1. The HIV proteins affect the Fas-mediated pathway in various ways. The picture illus- 
trates the type I and type II pathways. We consider the activation of Procaspase 3 to be the 
end of the apoptotic pathway (signifying the death of the cell). The type I pathway involves di- 
rect activation of Procaspase 3 by Caspase 8. The type II pathway requires signal amplification 
through the release of cytochrome c from the mitochondria resulting in Caspase 9 activating 
Procaspase 3. Tat upregulates Procaspase 8 and Bcl-2, but it can also down-regulate Bcl-2. Vpr 
upregulates Bcl-2 and downregulates Bax. HIV protease can cleave Bcl-2 into an inactive form 
and it can also cleave Procaspase 8 into active Caspase 8. 

5 Results of the Infected Model 

Based on the rules in Appendix A we simulated two different version of T cells. What we refer to 
as the "Non-latent" is a simulation of the T cell just after the contents of a virion have infected 
it, so there is Vpr active in the cytoplasm and the RNA is on its way to being incorporated into 
the DNA of the active CD4 cell. We call the second simulation method "Post-latent", which is 
a simulation of the reactivation of an infected, but latent CD4 cell. We seek to achieve greater 
understanding on the nature of latently infected cells. 

As we consider activation of Procaspase 3 to be an indication of the onset of apoptosis, we 
see in Fig. [2] that both the Post-latent and Non-latent cells die in approximately two days. Our 
simulations do indicate that reactivated latent cells will die more quickly than freshly infected 
cells. The lack of Vpr molecules in the latent cells is most likely the cause, as the anti-apoptotic 
effect of the protein would not be available. Also in Fig. [2] we show the truncation of Bid, a 
necessary step in the induction of the type II pathway. Active Caspase 8 is responsible for the 
truncation of Bid, so we again see the affects of the increased rate of Caspase 8 activation. The 
amplification of the signal through the mitochondria occurs earlier in the Post-latent than in 
the Non-latent simulation. 
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time (Hour) Time (Hour) 

Fig. 2. On the left, we see cell death for the latently infected cell occuring in 42 hours and 
freshly infected cell dying in 49 hours. On the right, we illustrate the rates of truncations for 
Bid, a signal for the induction of the type II pathway. 

In an HIV infected cell, Procaspase 8 can be activated by two ways, interactions with FADD 
or HIV Protease. The upregulation of Procaspase 8 by Tat increases the activation rate of 
Procaspase 8 (because it increases the chance of Procaspase 8 interacting with the existing 
FADD or HIV Protease molecules). As the number of molecules of Fas ligand is the same 
between both simulations, the increased rate of apoptosis is due to the creation and activity of 
HIV Protease and Tat. In Fig. [3l we show the activation of Procaspase 8 and Procaspase 9 for 
each of the simulations. 




time (Hour) Time (Hour) 

Fig. 3. These graphs show the increase of active Caspase 8 and Caspase 9 molecules over the 
entire simulation run. On the left is the Post-latent simulation and on the right is the Non-latent 
run. 

We next seek to determine whether the Procaspase 3 is being activated through the type 
I or the type II pathway. According to our rule set, an interaction between Procaspase 3 and 
active Caspase 8 or Caspase 9 can have two outcomes, resulting in activation of the Procaspase 
3 molecule or not. Interestingly, it appears that the two pathways work together to annihilate 
the cell. The Non-latent simulation shows the first interaction of a molecule of Procaspase 3 
occurs with a Caspase 8 molecule just after 10 hours into the run. It isn't until a 14 hours 
later (24.5 hours into the run, allowing for signal transduction through the mitochondria) that 
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we begin to see Procaspase 3 interactions with active Caspase 9. As we have shown in |21| . 
given a sufficiently high concentration of procaspase 8 in the T cell, there is no need for signal 
amplification. However, the initial level of Procaspase 8 in our model reflects Jurkat T Cells |19j . 
and we are seeing cooperation between Caspase 8 and Caspase 9 leading to cellular apoptosis. 
In Fig. H] we get an overall image of the effects of both pathways on Procaspase 3 activation. 




Time (Minutes) Time (Hour) 

Fig. 4. These graphs show the type I and type II pathways working together to activate Procas- 
pase 3 during the non-latent simulation. The Type I interactions are active Caspase 8 binding 
with Procaspase 3, and the Type II interactions are Caspase 9 binding with Procaspase 3. On 
the right we see the overall picture for the whole three days of simulation. On the left we can see 
a window of one whole minute from the simulation (from 26 hours to 26 hours and 1 minute). 

We find similarities in the Procaspase 3 activation between the Non-latent and Post-latent 
simulation runs. In the Post-latent run, we see type I interactions occur about 6 hours into 
the simulation, while type II molecular binding occurs after 17 hours (11 hours after the type 
I interactions). Fig. [5] show us an the overall interactions between Procaspase 3 with active 
Caspase 8 or Caspase 9. 




1560 _ 1561 6 12 18 24 30 36 42 48 54 60 66 72 

Time (Minutes) Time (Hour) 



Fig. 5. These graphs show the type I and type II pathways working together to activate Procas- 
pase 3 during the post-latent simulation. The Type I interactions are active Caspase 8 binding 
with Procaspase 3, and the Type II interactions are Caspase 9 binding with Procaspase 3. On 
the right we see the overall picture for the whole three days of simulation. On the left we can see 
a window of one whole minute from the simulation (from 26 hours to 26 hours and 1 minute). 
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The graphs in Fig. [4] and [5] imply that there are more interactions in the type I pathway. This 
does not necessarily mean that the type I pathway has a greater responsibility in activation of 
Caspase 3 than the type II pathway. Caspase 3 activation is a two step process. First, Caspase 8 
or 9 binds to the Procaspase 3 molecule, then after a period of time the binding ends with either 
an active Caspase 3 or the original Procaspase 3 molecule. Technically, based on the kinetic 
rates and the concentration of molecules involved, Caspase 8 can stay binded to Procaspase 3 
for longer than Caspase 9. Therefore, though it appears that Caspase 8 binds to Procaspase 3 
more frequently, it might not be the case that this results in more activation of Caspase 3 via 
Caspase 8 than Caspase 9. The graphs illustrate this possibility by the fact that you see constant 
jumps for the type II interactions (Caspase 9 and Procaspase 3), which implies less waiting time 
during the binding period. 

6 Discussion 

Based on the biological evidence from the literature, we were able to construct a simulation of 
the effects of HIV proteins on the Fas-mediated apoptosis pathway. We feel that our unique 
approach provides insight into the molecular dynamics of the Fas pathway. Our simulator leads 
to apoptosis in approximately two days, which is in agreement with observed results. But, we 
have shown that cooperation between the type I and type II pathway seems to be the method 
of the cell death. 

In this preliminary simulator, we have focused primarily on the difference between reactivated 
latent cells and actively infected cells, but we are looking to refine and expand our technique. 
We will apply new information on the nature of latent HIV infected cells as it becomes available. 
However we will also explore alternative HIV models. For instance, it would be interesting to 
model the effects HIV proteins on bystander cell apoptosis. The type 1 HIV has been implicated 
in reducing T cell levels by not only destroying the cell it invades, but through the killing of so- 
called bystander (uninfected) cells [14] . Various mechanisms have been debated as to the means 
through which the virus causes the destruction of the bystander (uninfected) cells. Along with 
Fas-induced apoptosis, other possible mechanisms for bystander cell death are reviewed in |31j . 
[29] and [2]. Upon being exocytised by an infected cell, several distinct types of HIV proteins 
can affect neighboring bystander cells. 

Soluble and membrane-bound Env can bind to the CD4 receptor of bystander cells (as is the 
case during infection). Along with subsequent binding to a co-receptor, it can induce apoptosis. 
In [10] and [5], it has been shown that ligation of the CD4 receptor by Env, is sufficient to 
increase apoptosis in bystander cells. The reasons for increased apoptosis following Env-CD4 
binding can be attributed to Bcl-2 down-regulation [18], increased Caspase 8 activation [lj, and 
upregulation of Fas [25] , FasL and Bax [31] . 

Extracellular Tat can be endocytised by a bystander cell, resulting in pro-apoptotic behavior. 
The addition of Tat to a culture of uninfected cells has been shown to increase apoptosis [23]. 
Endocytised Tat can up-regulate levels of Caspase 8 [4] and increase expression of the Fas ligand 
|31| . similar to its effects in infected cells. Extracellular Vpr can disrupt the mitochondrial 
membrane, leading to increased translocation of cytochrome c [31] . 

Lastly we would like to note the fact that the experimental information about the latent 
HIV infected T cells is scarce, the main reason being the fact that these cells are found in small 
numbers in vivo, thus it is rather hard to obtain sufficient numbers for experiments in vitro. We 
hope that our results will help guide the researchers in the lab and we are looking forward to 
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experimental results about these enigmatic cells that will suggest improvements and refinements 
of our model. 
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A The Model 



Our model consist of the following P system: 

n = (E, L, fi, M U M 2 ,M 3 , M 4 , M 5 , R u R 2 , R 3 , Ri, Rs) 

where: 

— In the alphabet we represent the proteins and complexes of proteins involved in the signalling 
transduction: 

S = {HIVRNA, Tat, Vpr, HIV Protease, FASL, FAS, FASC, FADD, FASC-FADD, FASC- 
FADD 2 , FASC-FADD3, FASC-FADD2-CASP8, FASC-FADD3-CASP8, FASC-FADD2-FLIP, 
FASC-FADD3-FLIP, FASC-FADD2-CASP82, FASC-FADD3-CASP82, FASC-FADD2-CASP8- 
FLIP, FASC-FADD3-CASP8-FLIP, FASC-FADD2-FLIP2, FASC-FADD3-FLIP2, FASC-FADD- 
CASP8, FASC-FADD-FLIP, CASP8, FLIP, FASC-FADD3-CASP83, FASC-FADD3-CASP82- 
FLIP, FASC-FADD 3 -CASP8-FLIP 2 , FASC-FADD3-FLIP3, CASP8f 41 , CASP8^, CASP3, 
CASP8^-CASP3, " CASP3*, 

CASP8^-Bid, tBid, Bid, Bax, tBid-Bax, tBid-Bax 2 , Smac, Smac*, Cyto.c, Cyto.c*, XIAP, 
Smac*-XIAP, Apaf, Cyto.c*-Apaf-ATP, CASP9, Cyto.c*-Apaf-ATP-CASP9, Cyto.c*-Apaf- 
ATP-CASP9 2 , CASP9*, CASP9*-CASP3, CASP9-XIAP, CASP3*-XIAP, Bcl2, Bcl2-Bax, 
Bcl2-tBid}l 

— The set of labels, L = {e, s, c, m, n}, can also be used to identify the different compartments 
where the signalling transduction takes place, namely, the environment (e), the cell surface 
(s), the cytoplasm (c), the mitochondria (to), and the nucleus (n). 

— In the membrane structure // = [1 [2 [3 [4 ]a[5 IsbbJij the compartments are represented. The 
compartments 1, 2, 3, and 4 can be labeled respectively with e, s, c, to, and n. 

— M{ = {wi,li} where the multiset Wj represents the initial amount of molecules in the com- 
partment and li is the label associated with the compartment. For Example: 

• h = e and M 1 = {FASL 301 } 

• l 2 = s and M 2 = {FAS 602 } 

• h = c and M 3 = {FADD 10038 , CASP8 20071 , FLIP 48779 , CASP3 120442 , Bid 15055 , 

Bax 50182 , XIAP 18066 , Apaf 60221 , ATP 6022141 CASP9 12044 , Vpr 700 } 

• / 4 = mandM 4 = { Smac 60221 , Cyto.c 60221 Bcl2 451661 } 

• h = n and M 5 = { HIVRNA 2 } 

— Ri, R2, R3, Ri, and R$ are the sets of rules associated with the environment, cell surface, 
cytoplasm, mitochondria, and the nucleus. These rules represent the chemical reactions that 
take place in each compartment or region of the cell. 

• Ri = {n} 

• R-2 = {r2,r4,r 6 ,r 8 ,r 1 o,rn,r 1 2,r li ,r 16 ,r 17 ,r 18 ,r 2 o,r22,r24,r26,r28,r30,r32,r3i, 

^36 , r 38 , r 40 , r 42 , r 44 , r 46 , r 48 , r 50 , r 52 , r 54 , r 56 , r 58 , r 60 , r 62 , r 63 , r 64 , r 65 , r 66 } 

• R3 = {^3,^5,r7,r9,rii,ri3,ri 5 ,ri 7 ,ri9,r2i,r23,r25,r27,r29,r3i,r33,r 35 ,r37,r39, 

r 4i , r 43 , r 45 , r 47 , r 49 , r 51 , r 53 , r 55 , r 57 , r 59 , r 6 i , r 67 , r 68 , r 69 , r 70 , r 71 , r 72 , r 73 , 
^74 , r 75 , r 76 , r 77 , r 78 , r 79 , r 80 , r 81 , r 82 , r 83 , r 84 , r 85 , r 86 , r 87 , r 88 , r 89 , r 90 , r 9 i , 
^92, ^93, r 94 , r 95 , r 96 r 9 8, r 9 9, rioo,noi, no2, no3, no4, Hos, »"106} 

• Ri = {r 97 } 

• -R5 = {no7,no8,no9,nio} 

2 Smac is an abbreviation for second mithocondria-derived activator of caspase. 
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We list the four sets of rules we simulated, which are based upon the apoptosis pathway. 
Each set incorporates a mechanism for Bcl2 involvement in the inhibition of apoptosis. 

label rule rate 

n : FASL[ FAS ] s -f [ FASC ] s k lf 

r 2 : [ FASC } s -> FASL[ FASC } s k lr 

r 3 : FASC[ FADD ] c -> FASC : FADD[ ] c k 2f 

r 4 : FASC : FADD[ ] c -> FADD ] c fc 2r 

r 5 : FASC : FADD[ FADD } c -> F^l^C : FADD 2 [ } c k 2f 

r 6 : F^IS'C : FADD 2 [ ] c -» F^IS'C : FAD£>[ FADD ] c fc 2r 

r 7 : F^,SC : FA£>£> 2 [ FADD ] c -> FAS'C : FADD 3 [ } c k 2f 

r 8 : F^45C : FADD 3 [ } c -» F^,SC : F^Z)Z) 2 [ F^DZ) ] c fe 2r 

r 9 : F^,SC : FADD 2 : C^,SP8[ F^DZ) ] c -» F^,SC : FADD 3 : C^,SP8[ ] c fc 2 / 

r 10 : F^SC : FADD 3 : C^,SF8[ ] c -> FAS'C : FADD 2 : CL4SP8[ FADD ] c fe 2r 

m : F^IS'C : FADD 2 : FLIP[ FADD ] c -» F^IS'C : FADD 3 : FLIP[ ] c k 2f 

n 2 : F^IS'C : FADD 3 : FLIP[ ] c -> F^S'C : FADD 2 : FLIP[ FADD ] c k 2r 

n 3 : F^5C : FADD 2 : C^5F8 2 [ FADD ] c -> FASC : FADD 3 : CA5F8 2 [ ] c fc 2/ 

n 4 : FA5C : FADD 3 : CA5F8 2 [ ] c -» FA5C : FADD 2 : CA,SF8 2 [ FADD } c k 2r 
n 5 : FAS'C : FADD 2 : CASP8 : FLIP[ FADD } c -> FAS'C : FADD 3 : CASP8 : FLIP[ } c k 2f 
n 6 : FAS'C : FADD 3 : CASP8 : FLIP[ ] c -» FAS'C : FADD 2 : CASP8 : FLIP[ FADD ] c /c 2r 

n 7 : FAS'C : FADD 2 : FLIP 2 [ FADD } c -» FAS'C : FADD 3 : FLFP 2 [ ] c fc 2/ 

n 8 : FASC : FADD 3 : FLIP 2 [ ] c -» FA5C : FADD 2 : FLIP 2 [ FADD } c k 2r 

n 9 : FA5C : FADD : CASP8[ FADD ] c -> FAS'C : FADD 2 : CASP8[ } c k 2f 

r 20 : FA5C : FADD 2 : CA,SF8[ ] c -> FAS'C : FADD : CASP8[ FADD } c k 2r 

r 21 : FASC : FADD : FLIP[ FADD ] c -> FAS"C : FADD 2 : FLIP[ } c k 2f 

r 22 : FASC : FADD 2 : FLIP[ ] c FASC : FADD : FLIP[ FADD ] c /c 2r 

r 23 : FASC : FADD 3 [ CASP8 ] c -> FASC : FADD 3 : CASF8[ ] c fc 2/ 

r 24 : FASC : FADD 3 : CASF8[ ] c -> FASC : FADD 3 [ CASP8 } c k 2r 

r 25 : FASC : FADD 3 [ FLIP } c -> FASC : FADD 3 : FLIP[ } c k 3f 

r 26 : FASC : FADD 3 : FLIP[ ] c -> FASC : FADD 3 [ FLIP } c k 3r 

r 27 : FASC : FADD 3 : CASF8[ CASF8 ] c -» FASC : FADD 3 : CASP8 2 { } c k 3f 

r 28 : FASC : FADD 3 : CASF8 2 [ ] c -» FASC : FADD 3 : CASP8[ CASP8 } c k 3r 

r 29 : FASC : FADD 3 : CASF8[ FLIP } c -» FASC : FADD 3 : CASF8 : FLIP[ } c k 3f 

r 30 : FASC : FADD 3 : CASP8 : FLIP[ ] c -» FASC : FADD 3 : CASP8[ FLIP } c k 3r 

r 31 : FASC : FADD 3 : FLIP[ CASP8 ] c -» FASC : FADD 3 : CASP8 : FLIP[ ] c fc 3/ 

r 32 : FASC : FADD 3 : CASF8 : FLIP[ ] c -» FASC : FADD 3 : FLIP[ CASP8 } c k 3r 

r 33 : FASC : FADD 3 : FLIP[ FLIP } c -» FASC : FADD 3 : FLIP 2 [ } c k 3f 

r 34 : FASC : FADD 3 : FLIP 2 [ ] c -» FASC : FADD 3 : FLIP[ FLIP } c k 3r 

r 35 : FASC : FADD 3 : CASF8 2 [ CASF8 ] c -> FASC : FADD 3 : CASF8 3 [ ] c fc 3/ 

r 36 : FASC : FADD 3 : CASF8 3 [ ] c -» FASC : FADD 3 : CASP8 2 [ CASP8 } c k 3r 

r 37 : FASC : FADD 3 : CASF8 2 [ FLIP ] c -> FASC : FADD 3 : CASF8 2 : FLIP[ ] c jfe 3/ 

r 38 : FASC : FADD 3 : CASF8 2 : FLIP[ } c -> FASC : FADD 3 : CASF8 2 [ FLIP } c k 3r 
r 39 : FASC : FADD 3 : CASF8 : FLIP[ CASP8 ] c -> FASC : FADD 3 : CASF8 2 : FLIP[ } c k 3f 
r 40 : FASC : FADD 3 : CASF8 2 : FLIP[ ] c -> FASC : FADD 3 : CASF8 : FLIP[ CASP8 } c k 3r 
r 4i : FASC : FADD 3 : CASF8 : FLIP[ FLIP ] c -> FASC : FADD 3 : CASP8 : FLIP 2 [ ] c jfe 3/ 
r 42 : FASC : FADD 3 : CASP8 : FL/F 2 [ ] c -» FASC : FADD 3 : CASP8 : FLIP[ FLIP ] c k 3r 

r 43 : FASC : FADD 3 : FL/F 2 [ CASF8 ] c -> FASC : FADD 3 : CASP8 : FLIP 2 [ ] c k 3f 

r 44 : FASC : FADD 3 : CASF8 : FL/F 2 [ ] c -» FASC : FADD 3 : FLIP 2 [ CASP8 } c k 3r 
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label rule rate 

r 45 : FASC : FADD 3 : FLIP 2 [FLIP ] c -> F^,SC : FADD 3 : FLLP 3 [ ] c fc 3/ 

r 46 : FASC : FADD 3 : FLIP 3 [ ] c -> FAS'C : FADD 3 : FLIP 2 [FLIP ] c k 3r 

r 47 : FAS'C : FADD 2 [ CASP8 } c -» F^5C : FAFF 2 : CA5P8[ ] c " % 

r 48 : FAS'C : FAFF 2 : C^15P8[ ] c -> FASC : FADD 2 [ CASP8 } c k 3r 

r 49 : FASC : FADD 2 [ FLIP } c -» F^5C : FADD 2 : FLIP[ ] c % 

r 50 : FAS"(7 : FADD 2 : FLIP[ ] c -» F^IS'C : FADD 2 [ FLIP ] c k 3r 

r 5 i : FAS'C : FAFF 2 : CA5P8[ ] c -> FAS'C : FADD 2 : C^5P8 2 [ ] c % 

r 52 : FAS'C : FAFF 2 : C^15P8 2 [ ] c -> FAS'C : FADD 2 : C^15P8[ CA5P8 ] c fc 3r 

r 53 : FAS'C : FAFF 2 : CASP8[FLIP ] c -> FAS'C : FADD 2 : CA5P8 : FLJP[ ] c % 

r 54 : FAS'C : FADD 2 : CAS PS : FLIP[ } c -> FAS'C : FADD 2 : CA5P8[ FLIP ] c fe 3r 

r 55 : PAS'C : FADD 2 : FLIP[ CASP8 ] c -> PAS'C : FADD 2 : (L4SP8 : PL/P[ ] c % 

r 56 : PAS'C : FADD 2 : CA5P8 : PL/P[ ] c -> PAS'C : : FLIP[ CASP8 ] c fe 3r 

r 57 : PAS'C : FADD 2 : FLIP[ FLIP ] c -> FASC : FADD 2 : FLIP 2 [ } c k 3f 

r 58 : FASC : FADD 2 : FLIP 2 [ } c -> FASC : FADD 2 : FLIP[ FLIP ] c k 3r 

r 59 : FASC : FADD[ CASP8 ] c -> FASC : FADD : CASP8[ ] c k 3f 

r 60 : FASC : FADD : CASP8[ ] c -> FASC : FADD[ CASP8 ] c fe 3r 

r 61 : FASC : FAFF[ FL/P ] c -> FASC : FADD : FLIP[ ] c k 3f 

r 62 : FASC : FAFF : FLIP[ ] c -> FASC : FAFF[ FL/P ] c /c 3r 

r 63 : FASC : FADD 2 : CASP8 2 { ] c -> FASC : FADD 2 [ CASP8^ 1 } c k 4 

r 64 : FASC : FADD 3 : CASP8 3 [ ] c -> FASC : FADD 3 : CASP8[ CASP8f 41 ] c fc 4 
r 65 : FASC : FADD 3 : CASP8 2 : FL/P[ ] c -» FASC : FADD 3 : FLFP[ CASP8f 41 ] c A; 4 

r 66 : FASC : FADD 3 : CASP8 2 [ ] c -> FASC : FADD 3 [ CASP8% 41 } c ' k 4 

r 67 : [ CASP8^ 41 } c -» [ CASP8^ ] c fc 5 

r 68 : [ CASP8 2 ,CASP3 ] c -» [ CASP8^ : CASP3 ] c fc 6/ 

r 69 : [ CASP8^ : CASP3 ] c -> [ CASP8^, CASP3 ] c fe 6r 

r 70 : [ CASP8^, CASP3* ] c -> [ CASP8^ : CASP3 ] c fc 7 

r 71 : [ CASP8^, Fid ] c -> [ CASP8^ : Fid ] c fc 8/ 

r 72 : [ CASP8^ : Fid ] c -» [ CASP8* 2 , Bid } c k 8r 

r 73 : [ CASP8* 2 , tBid } c -» [ CASP8^ : Fid ] c /c 7 

r 74 : [ tBid, Bax ] c — > [ tFiii : Fax ] c fcgj 

r 7 5 : [ tFiti : Fax ] c — > [ tBid, Bax ] c kg r 

rjQ : [ tBid : Bax, Bax ] c — > [ tBid : Bax 2 ] c kgj 

r 77 : [ tBid : Bax 2 ] c — > [ tFid : Fax, Fax ] c /cg r 

r 7 § : tBid : Bax 2 [ Smac] m — > Smac*[ ] m fcio 

r 79 : tFiti : Fax 2 [ Cyto.c ] m — > Cyto.c*[ \ m k\o 

r 80 : [ Smac*, X/AP ] c -» [ 5mac* : X/AP ] c jfe n/ 

r 84 : [ Smac* : XI AP } c -» [ Smac*, X/AP ] c fc llr 

r 82 : [ Cyto.c*, Apa/, AFP ] c -> [ Cyto.c* : Apaf : ATP ] c fc 12/ 

r 83 : [ Cyto.c* : Apaf : AFP ] c -> [ Cyto.c* , Apaf , ATP ] c k 12r 

r 84 : [ Cyto.c* : Apaf : ATP, CASP9 } c -» [ Cyto.c* : Apa/ : AFP : CASP9 ] c k 13f 

r 85 : [ Cyto.c* : Apaf : ATP : CASP9 ] c -> [ Cyto.c* : Apa/ : AFP, CASP9 ] c fci 3r 
r 86 : [ Cyto.c* : Apaf : AFP : CASP9, CASP9 ] c -> [ Cyto.c* : Apa/ : AFP : CASP9 2 ] c fci 4/ 
r 87 : [ Cyto.c* : Apaf : AFP : CASP9 2 ] c -> [ Cyto.c* : Apaf : ATP : CASP9, CASP9 ] c fc Ur 
r 88 : [ Cyto.c* : Apa/ : AFP : CASP9 2 ] c -»• [ Cyto.c* : Apa/ : AFP : CASP9, CASP9* ] c k 15 

r 89 : [ CASP9*, CASP3 } c -» [ CASP9* : CASP3 ] c fc 16/ 

r 90 : [ CASP9* : CASP3 ] c -> [ CASP9*, CASP3 ] c fci 6r 
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label rule rate 

r 91 : [ CASP9* : CASP3 ] c -> [ CASP9*, CASP3* ] c k 17 

r 92 : [ CASP9, XI AP } c -» [ CASP9 : IMP ] c A: 18/ 

r 93 : [ CASP9 : IJiP ] c -> [ CL45P9, XL4P ] c A: 18r 

r 94 : [ CASP3*,XIAP ] c -» [ C.4SP3* : IJiP ] c fcj 9/ 

r 95 : [ CASP3* : XL4P ] c -» [ CASPS*, IJiP ] c fei 9r 

r 9 6 : Pa:r[ PcZ2 ] m — ► [ PcZ2 : Bax ] m k 2 Qf 

r 9 7 : [ Bcl2 : Pax ] m — > Pax[ Pc/2 ] m k 2 o r 

r 9 g : tBid[ Pc/2 ] m — > [ Pc/2 : tBid ] m k 2 Qf 

r 99 : [ PcZ2 : tBid } m — > tBid[ Bcl2 ] m k 2 Q r 

Hoo : Vpr[ } m Vpr[ Bell } m k 21 f 

rioi : [ Vpr, Bax ] c -> [ Vpr ] c k 22f 

r W2 : [Tat] c ^ [Tat,Casp8] c k 23f 

rio3 : [ Tat } c -»• [ Tat, Bcl2 } c fc 24 / 

rio4 : [ Tat, Bcl2 } c -»■ [ Tai ] c fc 25 / 
rio5 : [ HIV Protease, Casp8 ] c — > [ HIV Protease, CaspS* ] m k 2 %f 

rio6 : HIVProtease[ Bcl2 ] m — > PP/Protease[ ] m fe 2 7/ 

r 107 : HIVRNA[ ) n -> [ HIVRNA } m k 28f 

r 108 : [ HIVRNA } n -► Fpr[ HIVRNA } m k 29f 

r 109 : PP/PAM[ ] n -► Tot[ HIVRNA } m k 30f 

mo : P/^PiVA[ ] n -> HIVProtease[ HIVRNA } m k 3 if 



The deterministic kinetic constants (reaction rates) mentioned in the previous table are 
given in the following; we refer the interested reader to [19] for more details about the rates and 
references for their estimation. 
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